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Abstract 



The equilibrium structure and exchange constants of CaMnOa and 
LaMnOs have been investigated using total energy unrestricted Hartree-Fock 
(UHF) and localised orbital configuration interaction (CI) calculations on the 
bulk compounds and Mn20j^~ and Mn20}^~ clusters. The predicted struc- 
ture and exchange constants for CaMnOa are in reasonable agreement with 
estimates based on its Neel temperature. A series of calculations on LaMnOs 
in the cubic perovskite structure shows that a Hamiltonian with independent 
orbital ordering and exchange terms accounts for the total energies of cubic 
LaMnOs with various spin and orbital orderings. Computed exchange con- 
stants depend on orbital ordering. Exchange contributions to the total energy 
vary between -20 and 20 meV per Mn ion, differences in orbital ordering en- 
ergy vary between 3 and 100 meV and a Jahn- Teller distortion results in an 
energy reduction of around 300 meV. The lattice constant of the lowest energy 
cubic perovskite structure (3.953 A) is in good agreement with the lattice con- 
stant of the high temperature 'cubic' phase of LaMnOs (3.947 A). The total 
energy of Pnma LaMnOs was minimised by varying lattice parameters and 
seven internal coordinates and a structure 194 meV per Mn ion below that 
of a structure determined by neutron diffraction was found. This optimised 
structure is nearly isoenergetic with a cubic perovskite structure with a 5 
percent Jahn- Teller distortion. UHF calculations tend to underestimate ex- 
change constants in LaMnOs, but have the correct sign when compared with 
values obtained by neutron scattering; exchange constants obtained from CI 
calculations are in good agreement with neutron scattering data provided 
the Madelung potential of the cluster is appropriate. Cluster CI calcula- 
tions reveal a strong dependence of exchange constants on Mn d e^ orbital 
populations in both compounds. CI wave functions are analysed in order to 
determine which exchange processes are important in exchange coupling in 
CaMnOs and LaMnOs- 
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I. INTRODUCTION 



CaMnOs and LaMnOs are endpoint compounds in the series Cai-^La^^MnOs, which 
has been thoroughly studied experimentally and theoreticallyS. They have relatively simple 
atomic and magnetic structures, their magnetic excitations are well described by a spin wave 
Hamiltonianlli and exchange constants, J, are well established by neutron scattering^'! and 
from the Neel temperatureafl. Exchange coupling in manganites has been extensively studied 
using model Hamiltonian&B and ab initio calculationsclii^. This paper presents results of 
bulk ab initio Unrestricted Hartree-Fock (UHF) and cluster Configuration Interaction (CI) 
calculations of exchange constants for both compounds. 

Exchange constants obtained from CI calculations are in excellent agreement with ex- 
periment and the localised orbital CI wavefunctions are analysed to determine which quan- 
tum fluctuations are most important in exchange coupling. Model Hamiltonian calcula- 
tions have attributed the exchange coupling, energy to to Mn^+Z^"*" superexchangei, 
Mn d^"*" d^^/Mn d^+ d^+ superexchang^, or bothl. Results of calculations presented 
below show that both exchange mechanisms operate and that O superexchange is the more 
important of the two. This was also found to be the case in the model Hamiltonian calcu- 
lations of Meskine et a/.i. 

CI cluster calculations provide detailed information on exchange couplings between neigh- 
bouring Mn ions, however a bulk electronic structure technique is required to study orbital 
ordering in LaMnOs. Total energies of LaMnOs with A and G-type anti-ferromagnetic 
(A-AF and G-AF) and ferromagnetic (FM) spin orderings have been computed in several 
isovolume structures in order to establish whether or not orbital ordering and spin ordering 
terms in the Hamiltonian for LaMnOs are independent. Obviously exchange constants will 
depend on orbital ordering, as the latter determines which empty orbitals are available to 
participate in exchange coupling. However, it is not known whether the e^, electron density 
in LaMnOs for a particular orbital ordering depends on spin ordering. It is shown below that 
a common orbital ordering energy for any of several orbital orderings can be identified and 
that this energy is independent of spin ordering, to a high degree. Spin and orbital ordering 
terms in the Hamiltonian are therefore independent, although orbital ordering determines 
the exchange constants. 

At low temperatures, CaMnOs exists in a cubic perovskite structure (lattice con- 
stant 3.73 A) with G-AF magnetic ordering^ and a Neel temperature of 130 K. Using 
the Rushbrooke-Wood formulae!, this Neel temperature implies an exchange constant J = 
6.6 meV. Note that throughout this work the spin Hamiltonian is of the form due to Domb 
and Sykes0 

H=EJ^.^ (1) 

<ij> ^ 

Si is a spin operator, S is the magnitude of the total spin for an ion and Jij is the exchange 
constant for a pair of ions. This form is adopted for the Hamiltonian as it is the same as 
that adopted in modelling spin wave dispersion in neutron scattering studiesi'i, except for 
a small Dzyaloshinsky-Moriya term. 
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At low temperatures the space group of LaMnOa is Pnmm. The ground state magnetic 
structure is A-AF and the unit cell contains four formula units consisting of rotated and 
distorted octahedra. There is one more d electron per Mn ion (c./. CaMnOa), which occupies 
an Gg orbital and induces a Jahn- Teller distortion in each MnOg octahedron, resulting in 
three distinct Mn-0 bond lengths of 1.91, 1.97 and 2.18 A0. The occupied e^ orbital is a 
linear combination of d^^-y^ and d^z^^r^ orbitals. The largest orbital component lies along 
the most elongated Mn-0 bond. The Pnma structure is shown in Fig. |I]. a, b and c axes 
referred to below are indicated on this diagram. The results of a number of neutron and 
x-ray scattering studies of the structure of LaMn03&^ over a range of temperatures are 
summarised in Ref.0. The Pnma structure can be viewed as containing planes of Mn^"*" ions, 
each joined to its in-plane neighbours by pairs of short (1.91 A) and long (2.18 A) Mn-0 
bonds. Each Mn^+ ion in a particular plane is coupled to Mn^"*" ions in planes immediately 
above and below by two Mn-0 bonds (1.97 A). The symmetry of the Pnma structure is 
such that there is one in-plane (nearest neighbour) exchange constant (Jy ) and one out-of- 
plane constant (J_l ). Both Jy and J_l have been determined from two independent neutron 
scattering studies to be -6.7 and 4.8 meV, respectively,ili. Thus there is FM coupling within 
planes and AF coupling between planes. 



II. EXCHANGE COUPLING MECHANISM 

The first comprehensive attempt to explain atomic and magnetic structures in doped 
and undoped manganites was made by Goodenough0 in 1955. He assumed three classes 
of exchange interaction between neighbouring Mn ions in undoped CaMnOs and LaMnOa 
lattices. A specific ordering of empty eg orbitals and relative orientations of pairs of empty eg 
orbitals correponding to Goodenough's classification are illustrated in Fig. ^ When empty 
eg orbitals are available on a pair of neighbouring Mn ions and are oriented towards one 
another (Fig. 0b) then AF coupling of the Mn ion spins is energetically favoured. This is 
because electrons from the central 0^~ ion of either spin are postulated to delocalise onto 
both Mn ions simultaneously, owing to the favourable exchange interaction (Hund's rule) 
between the delocalised electron and the Mn ion spin. However, if the Mn ion spins are FM 
aligned, only the electron from the central 0^~ ion with the same spin orientation as the Mn 
ions can delocalise onto either Mn ion, resulting in a higher energy for that state. Thus the 
empty orbital arrangement shown in Fig. results in an AF coupling of Mn spins. This 
is a type I exchange interaction according to GoodenoughB When one empty e^ orbital is 
suitably oriented for 0^~ ion electron delocalisation (Fig. |^), FM coupling of the Mn ion 
spins is favoured. This is a type II interaction. Finally, when no empty hybrids are available 
(Fig. ||d), no delocalisation occurs. This is a type HI interaction. This model has been used 
to explain the relative energies of A-AF, G-AF and FM magnetic states of CaMnOa and 
LaMnOa with a cubic perovskite structured. In that work it was found that the relative 
energies of these magnetic structures could be explained by counting the numbers of each 
type of interaction in each magnetic state and calculating the relative energy of each type 
of interaction. For both CaMnOa and LaMnOa it was found that the type I AF interaction 
was more energetically favourable than the type II interaction by 10 me The simplified 
description of exchange interactions just given assumes that an empty e^ orbital is either 
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available or not. However, empty Sg orbitals in LaMnOs are not purely d^^^yi or d^z^^r'^ in 
characteri. The mixed character of the empty e^ orbital permits some exchange coupling 
even when the empty e^ orbital is not optimally oriented. 

This type of reasoning was used by Millisi in a calculation of exchange coupling energies 
in CaMnOs and LaMnOs. In that work configurations allowed by the Pauli principle in 
which one or two electrons hop from the central ion to one or both neighbouring Mn 
ions are considered. Configurations which differ by a single electron hop have a single 
hopping matrix element, t. Diagonal elements of the Hamiltonian are parametrised using 
the energy required to excite one electron or a pair of electrons from an ion to a Mn 
ion. In CaMnOs the configuration which is assumed to lead to stabilisation of the AF state 
over the FM state is one in which a pair of electrons on the ion is excited onto separate 
Mn^+ ions. If this were indeed the origin of exchange coupling in CaMnOa then one would 
expect this configuration to appear in an ab initio ground state CI wave function, but this 
is not the case. However, the main idea of this model, that more low energy configurations 
are available to singlet states than high spin multiplicity states, is in accord with results 
presented here. 

CI cluster calculations of exchange constants in La2Cu04ii@ and KNiFail, which used 
delocalised molecular orbitals have been reported quite recently. CI cluster calculations 
described below were carried out in a localised orbital basis. The localised orbital basis 
provides a means of identifying the exchange coupling mechanism in terms of fluctuations 
of electrons between localised orbitals. These calculations were performed on Mn20ii~ and 
Mn20j5~ clusters representing fragments of CaMnOa and LaMnOs. Details of the calcu- 
lations, including the method used to generate the localised orbitals, details of a spherical 
array of point charges surrounding the clusters, etc. are given in Appendix 0. The wave- 
functions for the clusters contain orbitals which are partitioned into a (doubly-occupied) core 
orbital space, an active space containing the 2p orbitals of the ion situated between the 
two Mn ions in the cluster as well as singly-occupied Mn d orbitals and an external space 
containing unoccupied Mn d orbitals. The core orbital space contains 'core' electrons as well 
as valence electrons not in the active orbital space. The ions in the clusters treated quantum 
mechanically consisted of two corner sharing MnOe octahedra. The localised orbitals in the 
active space for the Mn20i^~ and Mn20if~ clusters are shown in Figs. ^ and respec- 
tively. The main electronic configuration for the Mn20jf~ cluster representing CaMnOa is 
one in which each Mn ion with a formal 4-f charge contains three t2g electrons and each 
oxygen ion exists in a closed shell configuration. The actual charge on the Mn ions is 
significantly reduced as there is a covalent component to the Mn-0 bonds, as can clearly 
be seen in the contour plot of the localised orbital with mainly O 2pz character in the top 
panel of Fig. ^. The actual Mn ion charge in CaMnOs is +2.13, according to a MuUiken 
population analysis of the UHF wave functions obtained for CaMnOs. The formal charge on 
Mn ions in LaMnOs is 3+, however a MuUiken population analysis of UHF wave functions 
for LaMnOs yields a charge of +2.24 The O ion charges in the two compounds are -1.33 
(CaMnOs) and -1.75 and -1.82 (LaMnOs). Exchange constants were calculated by finding 
the energy difference between the spin singlet and spin septet(nonet) states of the Mn20ii~ 
and Mn20}i~ clusters. 

Wave functions were constructed from the localised orbitals shown in Figs. |^ and ^ and 
doubly occupied core orbitals. A septet state for the Mia20\\~ cluster was constructed from 
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six singly-occupied t2g orbitals and doubly occupied O 2p orbitals localised on the central 
O ion in the cluster. The form of this wave function is 



^septet ^ A({core}(0a;y,i0x^,/</>y^,«0a'j/,r0a'^,r0?;^,r)(aaaaaa)) (2) 

A is the anti-symmetrising operator and the subscripts / or r on t2g orbitals in Eq. |^ 
indicate that they are centred on the left or right Mn ion, respectively, {core} is a product 
of doubly occupied orbitals in the core orbital space which includes the 2p orbitals on 
the central 0^~. This is the restricted open shell Hartree-Fock (ROHF) wave function for 
the septet state, constructed using localised molecular orbitals. Self-consistent field (SCF) 
ROHF wave functions can be computed using a number of standard electronic structure 
packages such as the GAMES S packageEa used in this work. 

The singlet state is constructed from the same set of singly-occupied orbitals with a 
spin-coupling of the form 

-j={aaappp — PPpaaa) (3) 
v2 

This is one of five spin eigenfunctions^ (SEE) for six electrons coupled into a singlet 
state. Provided that the spatial orbitals multiplying this SEE are ordered such that orbitals 
localised on each Mn ion are grouped together, we expect this SEE to dominate the CI wave 
function, since Hund's rule requires spins on each ion to be coupled with the same spin. 
This is indeed found to be the case in the actual CI wave function for the singlet state of the 
Mn20jf~ cluster representing CaMnOs. The wave function for the singlet state is therefore 

^singlet ^ -^A(j^core}{(j)^y^i(j)^,^i(j)y,^i(j)^y^r(pxz,Ayz,r^ {aaa(3f5(5 - f3P(3aaafj (4) 

Using conventional rules for evaluating determinantal energiesii, the energy difference 
between the two states is Kxz,i xz,r + ^yz,i yz,r, with the singlet state lying above the septet 
state (assuming that other inter-site exchange integrals are zero because of negligible spatial 
overlap). When the ground state energies of the singlet and septet states of the Mn20ii~ 
cluster with wave functions in Eq. |^ and ^ were evaluated, the singlet state was 3.6 meV 
above the septet state. This implies a value of 1.8 meV for the exchange integrals just men- 
tioned. Note that we use the notation Kjj for exchange integrals between specific molecular 
orbitals while we use the notation J for the (effective) exchange coupling energy of two spins 
on different Mn ions. The singlet and septet states of this configuration are analogous to 
the Heitler-London valence bond wave function for the singlet and triplet states of the He 
atom in a ls2s configuration. In that case the triplet state is lower than the singlet state by 

In general, CI wave functions with N electrons in the active orbital space consist of linear 
combinations of spin-adapted functions (SAF's) 

^c7 = Ec.^P^ (5) 

i 

ijf^^ = A({core}(^>fc...0,0ie,) (6) 
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where a SAF is a product of spatial orbitals and a SEF, 9a, for the particular spin state 
in question. The septet and singlet SAP's in Eq. § and ^ are the dominant terms in a more 
general CI expansion of the septet and singlet wave function of the Mn20i^~ cluster. All 
SAF's which are obtainable by exciting one or two electrons from the dominant SAP's to 
empty orbitals in the active space are included in the expansion. As stated above, the 13 
orbitals in the active space in the calculations described here are comprised of 10 orbitals of 
mainly Mn 3d character and 3 of mainly O 2p character localised on the O ion between the 
two Mn ions. These excited electron SAP's enter the wave function with a maximum weight 
of order 10~^ and a corresponding occupancy of order 10~^ and it is these which lower the 
energy of the singlet state below the septet state when the spins are AP coupled. The main 
excited SAP's in the singlet and septet wave functions are those in which: one electron is 
excited from an O 2p orbital to the Mn e^ orbital aligned with the Mn-0 axis (O to e^ (le)); 
a pair of electrons are excited from one O 2p orbital to the same Mn e^ orbital (O to e^ 
(2e)) and an excitation in which a t2g electron is transferred fron one Mn ion to the other 
{t2g exchange). Obviously the latter exchange process is only allowed in the singlet state 
as it violates the Pauli exclusion principle in the septet state when the t2g shells are half 
filled, as in CaMnOa. Excitations in which a pair of electrons are excited from the O ion to 
separate Mn ions are found to have negligible weight for both spin states. 



III. RESULTS 

A. CaMnOs: Bulk UHF Calculations 

UHP total energy calculations were performed using the CRYSTAL98 co d(S for PM, A-, 
C-, and G-AP spin orderings. The energy of the cubic PM structure with the experimental 
lattice constant of 3.73 A was adopted as the reference energy (0 meV); calculations were 
also performed for each of the spin orderings with a lattice constant of 3.75 A. Total energies 
and magnetic moments from these calculations are given in Table |. When these total energy 
differences are fitted to the Hamiltonian in Eq. |l| with nearest (Ji) and second nearest (J2) 
neighbour interactions {i.e. along [a, 0,0] and [a,a,0], where a is the lattice constant), the 
parameters obtained for a lattice constant of 3.73 A are Ji = 10.7 meV, J2 = 0.3 meV and 
for a lattice constant of 3.75 A, the parameters are Ji = 10.1 meV, J2 = 0.3 meV. It is 
generally believed that exchange interactions which connect magnetic ions along a linear 
chain are stronger than those which do not, such as the J 2 interaction here. However, in 
the cubic perovskite structure, exchange interactions along [2a,0,0], etc. contribute equally 
to all four spin orderings studied and so cannot be extracted from the data presented here. 
Similar values for Ji have been obtained from model Hamiltonian calculations by Meskine 
et al. (Ji = 9.5 meV)i. Note that the definition used for the exchange energy in that work, 
the difference between the energy of a pair of ferro- and anti-ferromagnetically coupled Mn 
ions, is twice the exchange energy defined in Eq. |l] above. Hence values for exchange energies 
from that work have been divided by two in order to compare them to values in the present 
work. 
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B. CaMnOs: Cluster CI Calculations 



Exchange energies obtained from cluster CI calculations depend strongly on Mn e^ and 
O 2p orbital populations. In turn these populations depend on the Madelung potential of 
a sphere of point charges surrounding the Mn20jf~ cluster. The charges were located on 
crystal ion sites and MuUiken populations of ions in bulk UHF calculations on CaMnOs were 
used as a guide in choosing the magnitudes of these charges. The sphere of point charges 
had a radius of just over 20 A and contained around 3300 charges. The radius was chosen 
so that the sphere was overall almost charge neutral; each unit cell of point charges was also 
neutral. The dependence of exchange energies in CaMnOa and LaMnOa on Mn ion charge, 
measured by the Mulliken population of that ion in the SCF cluster calculation, is shown in 
Fig. ^ The charge on the two Mn ions in the cluster was adjusted by transferring charge 
from Mn point charge sites to Ca or La point charge sites in the sphere of point charges. 
The total charge of the Mn and Ca (or La) point charges was kept constant and the O 
ion charge was maintained at the UHF Mulliken population value. It can be seen that the 
magnitude of the exchange energy increases as charge is removed from the Mn ion, which is 
reasonable as charge is mainly being transferred to/from the e^ orbitals which are directly 
involved in the exchange coupling mechanism - as the e^ orbital becomes filled, the exchange 
energy diminishes. The CI calculation value of Ji = 8.1 meV quoted for CaMnOa in Table 
|I|is the value obtained for an Mn cluster ion charge of +2.13, the Mn ion charge determined 
from the UHF calculation. This is to be compared to an estimate of the experimental value 
of Ji = 6.6 meV, derived from the Neel temperature of CaMnOs. 

The fundamental SAF's for the septet and singlet states of the Mn20ji~ cluster were 
given in Eq. ^ and |^. In the fundamental SAF wave functions for either spin state, each has 
a SAF coefficient, Cj, of unity, however when additional SAF's are permitted in the wave 
function {i.e. permitting ion 2p to e^ excitations, etc.) the weights of fundamental 
SAF's are around 0.9950 and additional SAF's corresponding to O superexchange and i2g 
exchange enter the wave function with SAF coefficients of order 0.01. Even for limited active 
spaces (as in these calculations) the number of SAF's entering the wavefunction means that 
a convenient way to analyse the wavefunction is to tabulate the summed occupancies {i.e. 
\ cj I) of configurations of a particular type. There are, for exarnple, several SAF's in which 
one electron is excited from an O 2p orbital to an Mn e^ orbitalEJ. The relative magnitudes 
of these occupancies are a measure of the importance of each type of fluctuation about the 
fundamental SAF conflgurations. Summed occupation numbers for the Mn20i^~ cluster are 
given in Table |ITT[ It can be seen that the fundamental (or main) SAF has occupancy 0.9926 
for the singlet state while it has occupancy 0.9943 in the septet state and therefore that there 
are larger correlation effects in the singlet state. SAF's in which a electron has hopped 
from one Mn ion to the other have an occupancy of 0.0005, while these fluctuations are absent 
from the septet state owing to the Pauli exclusion principle, as noted above. However, the 
main difference in septet and singlet wave functions is in the occupancy of states in which 
one electron is transferred from an O 2p orbital to an e^ orbital, the occupancy being 0.0038 
for the singlet state and 0.0027 for the septet state. The occupancy of SAF's in which a 
pair of electrons is transferred from O 2p to one Mn e^ orbital is the same for both spin 
states. The energies of both spin states relative to the energy of the fundamental septet 
SAF are also given in Table |T|. The septet state with O superexchange fluctuations is 133.4 
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meV below the fundamental septet SAF. This is the correlation energy for that statecj. The 
singlet state with O superexchange and t2g fluctuations is 149.6 meV below the reference 
energy and 153.2 meV below the fundamental singlet SAF energy. The latter energy is the 
correlation energy for the singlet state. Correlation energies for the Mn20i^~ and Mn20}f~ 
cluster CI wave functions are illustrated schematically in Figure ^. Correlation energies are 
around 50 percent larger in Mn20}^~ than in M'd.201i~ and this is reflected in the larger 
exchange energy in CaMnOs. It is worth noting that when the CI cluster calculation for 
the exchange energy in CaMnOa was performed with no point charge array surrounding 
the cluster, the exchange energy obtained was 57 meV, well in excess of the experimental 
value. This emphasises the importance of Madelung terms in the crystal Hamiltonian in 
determining exchange energies in strongly correlated materials. 



C. LaMnOs: Bulk UHF Calculations 

Total energy calculations were performed on LaMnOs in the ideal perovskite (cubic) 
structure, a tetragonal perovskite structure, a cubic structure with a Jahn- Teller distor- 
tion of the MnOg octahedra and the Pnma structure with atomic coordinates derived from 
experimentEl and by minimising the total energy by varying lattice parameters and internal 
coordinates not determined by symmetry. These structures are summarised in Table 0. 
The Jahn- Teller distortion consisted of elongation or contraction of Mn-0 bonds parallel to 
the ac axes of the unit cell. These are the Mn-0 bonds which induce FM coupling between 
Mn ions in the Pnma structure. The cubic structure with the lowest energy had a lattice 
constant of 3.953 A (volume 61.77 A'^ per Mn ion) which is comparable to the lattice constant 
of the 'cubic' phase of LaMnOs (3.947 A) which occurs at temperatures above 75OK0. All 
relative energies and lattice volumes will be assumed to be per Mn ion hereafter. When this 
structure is changed by a 5 percent Jahn- Teller distortion (Table ^V]) the energy is lowered 
by 304 meV and the magnetic ground state of the structure switches from (ix-2-y2 d^z^-r'^ 
FM to dsa;2-r2 dsj/2_r2 FM (sec below). 

The total energy of the Pnma structure using coordinates from experiment^ (Table is 
200 meV above the Jahn- Teller distorted structure. The total energy of the Pnma structure 
was mmimise by varying the lattice parameters and 7 internal coordinates not determined 
by symmetry of the Pnma space group. The total energy of the energy minimised structure 
is 6 meV above the Jahn- Teller distorted structure. The optimised lattice parameters and 
internal coordinates are given in Table |rV|; the a lattice vector is essentially unchanged while 
the b and c lattice vectors increase in magnitude by 1.1 and 1.6 percent, respectively. The 
lattice volume rises from 60.89 to 62.53 A^. Probably the most important changes which 
occur on minimising the total energy are: the degree of Jahn- Teller distortion is reduced; 
La-0 distances increase significantly. In the lowest energy cubic structure there is one Mn- 
O distance of 1.976 A and an La-0 distance of 2.795 A. On introducing the 5 percent 
Jahn- Teller distortion these become Mn-0 distances of 1.877, 1.976 and 2.075 A and La-0 
di.tance= ot 2.795 and 2.797 A. In the experimental Pnma *uctnre0 the Mn-0 distance, 
are 1.903, 1.957 and 2.185 A and the La-0 distances are 2.433, 2.461 and 2.548 A. These 
change to 1.910, 1.944 and 2.135 A and 2.609, 2.666 and 2.684 A in the energy minimised 
structure. Hence lower energies are found for structures with larger La-0 distances and a 
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reduced Jahn- Teller distortion. The combined ionic radii of La^"*" and are 2.76 Aca. La- 
O distances in the energy minimised Pnma and cubic structures lie just below the combined 
ionic radii distance, whereas the La-0 distances in the experimental Pnma structure lie 
well below this distance. The cubic structure with a Jahn- Teller distortion and the energy 
minimised Pnma structure are both lower in energy than the lowest energy cubic structure 
by around 300 meV. This energy lowering by a Jahn- Teller distortion is half of the lowering 
assumed by MiUiH in a calculation of electron-phonon coupling in Cau^^La^^MnOs. The 
UHF calculations reported here are similar to those reported by Su et a/.t3. They report an 
energy lowering of 1055 meV when the cubic structure is changed to the experimental Pnma 
structure with no volume change. This calculation will overestimate the energy difference 
between such structures as the cubic structure with the Pnma structure equilibrium volume 
is not the mimimum energy cubic structure. 

For the cubic perovskite structure it was found that variations of the total energies of 
different spin and orbital orderings can be fitted very well by a Hamiltonian of the form 

H = J^^^^ + Hoc (7) 

where Hqo is an orbital ordering term which depends only on the orbital order. For these 
calculations the cubic unit cell was doubled along [110], [101] and [Oil] directions (G-AF 
spin and orbital ordering) and along the [001] direction (A-AF spin and orbital ordering) 
and total energies and charge density difference plots§ were computed for dx2_y2 d^2_„2 , 
da.2_y2 d322_^2 and d3z2_,^2 d^^^-r'^ orbital orderings and FM, A-AF and G-AF spin orderings. 
The d^.2_^2 d322_r2 A-AF combination is incompatible with the unit cell doublings chosen 
and was omitted. Total energies are given in Table and charge density difference plots for 
each of these orbital orderings are shown in Fig. ^ 

For dx2-y2 dx2-y2 and d322_^2 d^z^^r^ orbital ordering, distinct exchange constants in the 
xy plane, Jy, and in the xz plane, J_l, are postulated, whereas for dj.2_y2 d322_^2 ordering a 
single exchange constant, J = Jy = J^, is postulated. Exchange constants for each orbital 
ordering are given in Table |V|. AF exchange constants are obtained when adjacent Mn 
orbital ordering is the same and FM coupling is observed when adjacent Mn e^ orbitals 
differ. This observation also applies to Pnma structures studied: FM coupling is observed 
between in-plane Mn ions with alternating e^ orbital orientations; AF coupling is observed 
when adjacent Mn e^ orbitals have the same orientation. Magnitudes of AF couplings vary 
between 5.1 and 14.2 meV and one FM coupling of -6.0 meV is observed in the dx^-y^ d^z^^r^ 
A-AF ordering. 

Once exchange constants have been computed for a particular orbital ordering, compar- 
ison of structures with the same magnetic structure but different orbital ordering permits 
differences in orbital ordering energies to be calculated. The actual magnitude of orbital 
ordering energy, Hqo, of course depends on the reference energy chosen. The choice of the 
FM dx2-y2 dx2-y2 structure as the reference energy structure yields values of -17.4, -20.0 and 
-113.4 meV for the da;2_2^2 dx2-y2 , d^z^-r^ d^z^^r^ and dx2_y2 d322_^2 relative orbital ordering 
energies. The important result here is that an alternating orbital order {dx2-y2 d^z^-r^ ) 
is around 90 meV below those with the same orbital order on each site (da.2_y2 d^.2_j^2 or 
d^z'^-r'^ d3z^_r2 ) in the cubic perovskite structure. When the values of Hqo and exchange 
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constants just mentioned are used to compute the relative energies of the eight spin and 
orbital orderings considered, the maximum deviation from the relative energies reported in 
Table |V| is 0.2 meV, demonstrating the suitability of the Hamiltonian in Eq. |^. The fact 
that charge density difference plots for different spin ordering and the same orbital ordering 
are very similar suggests that this should be the case. 

Using the fact that orbital and spin contributions to the Hamiltonian are independent, 
differences in total energy of a particular spin order as a function of lattice distortion may be 
attributed to differences in orbital ordering energy. Fig. ^ is a plot of total energy for each 
orbital ordering with G-AF magnetic order as a function of isovolume, tetragonal lattice 
distortion. These calculations were performed using P4/mmm space group symmetry. It 
can be seen that d^i-y^ d^z^^r^ orbital ordering is the most stable ordering only within a 
small parameter range about the cubic structure. When the tetragonal distortion is such 
that the lattice is elongated along the z axis, d^z^^r^ d^z^-r^ ordering is favoured and when 
it is compressed along this axis, dj.2_j^2 d^.2_j^2 ordering is favoured. This may be explained 
by a simple electrostatic argument - the ordering which is favoured in either case is the 
one where the occupied eg orbitals are oriented along the elongated axis or axes, thereby 
reducing the Coulombic repulsion energy. The greatest stabilisation relative to the cubic 
lattice is found for an x/z ratio of 0.94 where the energy is 164 meV below that of the cubic 
G-AF reference energy. This stabilisation is still significantly less than the stabilisation of 
298 meV which results when the energy minimised Pnma structure is adopted. 

Relative energies and exchange constants for the Jahn- Teller distorted structure and both 
Pnma structures studied are given in Tables [VII] and [VIII| , respectively. Charge density 
difference plots for the Jahn- Teller distorted and energy minimised Pnma structures are 
shown in Fig. ^. The magnetic ground state of the Jahn- Teller distorted structure is FM, but 
this is almost isoenergetic with the A-AF structure. This is because the in-plane exchange 
constant is FM while the out-of-plane exchange constant is FM (but small). The magnetic 
ground state of the cubic structure is G-AF with AF coupling between all neighbouring 
Mn ions. The switch to FM coupling between neighbouring Mn ions in-plane is due to the 
Jahn- Teller distortion in-plane. Both Pnma structures studied have A-AF magnetic ground 
states (as is the case in nature) but magnitudes of exchange constants obtained from these 
calculations are smaller than those obtained from neutron scattering datai'i (Table [VIII] ). 
Values of 0.6 and -3.7 meV for J^ and Jy may be compared to 0.8 and -3.5 meV obtained in 
a similar UHF calculationlli and 4.8 and -6.7 meV from experimentiH. A local spin density 
approximation (LSDA) calculation^ found values of -3.1 and -9.1 meV for J^ and Jy . This 
calculation did find an A-AF ground state for Pnma LaMnOa, however, as second nearest 
neighbour exchange constants are significant in the LSDA calculation and favour an A-AF 
magnetic ground state. 



D. LaMnOs: Cluster CI Calculations 

Cluster CI calculations for LaMnOs were performed using Mn20ii~ clusters with the 
Mn ions in the same configuration as a pair of Mn ions in the ac plane (Fig. |^) and with 
the Mn ions along a line parallel to the b axis. The former cluster corresponds to a pair 
of Mn ions which is expected to be FM coupled while the latter corresponds to a pair of 
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ions which is expected to be AF coupled. Mn e^ orbital ordering in the former cluster had 
the form illustrated schematically in Fig. ^ while the latter had orbital ordering as in Fig 
Clusters and surrounding point charges with the experimental Pnma structureEl and 
the energy minimised structure were used. Exchange constants for LaMnOs derived from 
these cluster calculations are given in Table |V111| . Cluster CI calculations with Mn, O and 
La surrounding point charges of 2.6, -1.8 and 2.8 (close to Mulliken population values from 
UHF calculations) result in exchange constants of 3.3 and -3.6 meV for and Jy when the 
experimental structure is used. These values change to 5.1 and -7.4 meV when the energy 
minimised structure (Table |Vp is used. 

The Madelung potential has an important role in determining exchange constants in 
manganites. Obviously ions several lattice constants or more distant from the ions in the 
central cluster may be treated as point charges rather than distributed charges without 
significantly altering the potential within the central cluster. However point charges adjacent 
to the central cluster may cause a significantly different potential within the cluster and affect 
the results of the exchange constant calculation. This question has previously been addressed 
by other workerJ^llii. In order to estimate the effect of terminating the cluster with point 
charges, cluster CI calculations were performed with the 12 La point charges nearest to the 
central cluster ions replaced by La^"*" pseudopotential£il. This resulted in a small increase 
in and no change in Jy compared to the calculation where only point charges were used. 
The values obtained for and Jy from these calculations were 5.2 and -7.4 meV, which are 
in good agreement with the experimental values: 4.8 and -6.7 meV. Values for the exchange 
constants derived from the model Hamiltonian calculations of Meskine et a/.i are also given 
in Table [VITT| . 

Relative energies and SAF occupancies for the Mn20j5~ clusters used for calculating 
exchange constants in LaMnOa in the energy minimised structure are given in Table ffX. 
The fundamental SAF singlet states are 11.9 meV (J_l) and 17.9 meV (Jy) above the nonet 
states of the clusters. When additional SAF's are permitted in the wave function the sin- 
glet(nonet) states are lowered by 105.4(83.3) meV (J_l) and 82.8(79.7) meV (Jy). These are 
the correlation energies for these states. The singlet state of the cluster used to calculate 

is 10.2 meV lower in energy than the nonet state giving a value for of 5.1 meV while 
the nonet state of the cluster used to calculate Jy is 14.7 meV lower in energy than the 
singlet state giving a value of -7.4 meV for Jy. From Table ^ it can be seen that O to e^ 
(le) excitations are the main fluctuations about the fundamental SAF state. The weight of 
the fundamental SAF in the singlet states of either cluster is less than in the nonet states, 
reflecting the greater degree of correlation in the singlet states. In the calculation the 
greater correlation energy of the singlet state c.f. the nonet state is sufficient to make the 
singlet state the ground state and give an AF exchange constant. On the other hand, in 
the J|| calculation the singlet correlation energy is just greater than that of the nonet state 
and, together with the fact that the singlet state of the fundamental SAF wave function is 
17.9 meV above the nonet state, this results in a nonet ground state and an FM exchange 
constant. 
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IV. DISCUSSION 



UHF and CI cluster calculations for the exchange constant in CaMnOs are in reasonable 
agreement with estimates for its value based orr the Rushbrooke-Wood formula! a,rd its 

Neel temperature. The calculated exchange constants are larger than the estimate based on 
experiment. The single AF exchange constant is mainly a result of O to e^ (le) excitations 
which lower the energy of the singlet state of a pair of adjacent Mn ions below that of the 
septet state. The magnitude of the exchange constant derived from CI cluster calculations 
depends strongly on the Madelung potential within the cluster and there is agreement be- 
tween theory and estimates based on experiment only when that potential results in ionic 
charges in the cluster similar to those in the bulk UHF calculation. 

LaMnOs is more complex than CaMnOs. It is also more ionic than CaMnOs with 
Mulliken populations of ions closer to the formal ion charges. A number of orbital and spin 
ordered states exist within a small energy range, say 300 meV, close to the ground state. 
The energies of several spin and orbital ordered states of cubic LaMnOa are well described 
by the Hamiltonian in Eq. 0. In the remainder of this section exchange constants in cubic 
and Jahn- Teller distorted LaMnOs are correlated with Mn^+ ion orbital ordering and 
ion charge density distortions and Mn^"*" ion interactions are identified as type I, II or III 
according to Goodenough's schemed. Finally, the role of correlation and availability of 
empty orbitals on magnetic ion sites in AF and FM coupling is discussed. 

Cubic LaMnOa has a dx'^-y'^ ^zz'^-r'^ FM ground state and a lattice constant of 3.953 A. 
Exchange constants depend on orbital ordering and range from -6.0 to 14.2 meV. Charge 
density difference plots (Fig. |^) show that the charge density on an Mn ion is essentially 
independent of charge densities on neighbouring ions. That density is determined solely 
by the ion's orbital ordering. However, charge densities on more polarisable ion sites 
depend on charge densities at both neighbouring Mn ion sites. For (i^z2_^2 (\^^2_^2 orbital 
ordering (Fig. 0, top panels), ions in Mn-0 bonds in the xy plane undergo a quadrupolar 
distortion in which charge is displaced from the Mn-0 bond axis into directions perpendicular 
to the bond while ions in Mn-0 bonds along the z axis are much less severly distorted 
and the ions tend to elongate along the bond axes. The (l^z2_j.2 character of the ordered Mn 
Gg orbitals can be seen clearly in the top right panel of Fig. 0. The exchange constant for 
Mn^+ ions in the xy plane with this orbital ordering is Jy = -0.1 meV while the exchange 
constant for Mn^+ ions along the z axis is 14.2 meV. Thus a weak exchange coupling is 
associated with the quadrupolar distortion of charge away from the bond axis while a much 
stronger coupling is associated with a nearly spherical ion in which charge density tends to 
concentrate along the bond axis, compared to the spherical 0^~ ion. 

For dx2_y2 dj.2_y2 orbital ordering the d^2_j^2 character of the ordered Mn e^ orbitals 
is clearly seen in the middle left panel of Fig. |^. There is a relatively weak quadrupolar 
distortion of the 0^~ ions in the xy plane and a stronger quadrupolar distortion of 0^~ ions 
along the z axis. The exchange constant for Mn^"*" ions in the xy plane is 5.1 meV while it 
is 7.2 meV for Mn^"*" ions along the z axis. 

The dx2_y2 or d3z^_r2 character of orbital ordering can be seen in both bottom panels 
in Fig. 0. In the xy and xz planes the 0^~ ion charge density is polarised in a circulating 
pattern, even though the 0^~ ions are situated midway between the Mn^+ ions. Charge 
is polarised towards regions at Mn3+ ion sites where there is a reduction in charge density 
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below that of spherical Mn^"*" ions, as indicated by negative contours in the charge density 
difference plots. Around each Mn^"*" ion with the 6.3^2 -r^ eg orbital occupied, charge is 
deformed towards the d^z^-r^ ion in the xy plane and away from it along the z axis, whereas 
for Mn^+ ions with the da.2_j^2 e^ orbital occupied, charge is deformed towards it along the z 
axis and away from it in the xy plane. Thus each Mn^^ ion is coupled to each neighbouring 
ion by a polarised O ion and there is one FM exchange constant of -6.0 meV. 

A simple pattern of orbital ordering is obtained for the Jahn- Teller distorted structure 
where the unit cell was doubled in the [001] direction (Fig. This pattern of orbital 
ordering was obtained without biasing the initial guess wave function in anyway (Appendix 
A). Orbital ordering in the xy plane is an alternating 6.3^2 d3j,2_^2 pattern which is repeated 
with period one along the z axis. This is the a-type orbital ordering discussed in Ref.il. 
The FM exchange constant between Mn^+ ions in the xy plane is -8.1 meV and the weak 
FM exchange constant along the z axis is -0.1 meV. There is strong deformation of 
ion charge density in the xy plane towards regions of reduced charge density at the Mn^+ 
ion sites, associated with strong FM exchange coupling (Fig. ^, left panel). There is strong 
quadrupolar distortion of the charge density at ion sites coupling Mn^"*" ions along the 
z axis, associated with a weak FM exchange coupling. 

The same a-type orbital ordering is also found in the Pnma structures studied. There is 
a similar pattern of circulating charge polarisation towards regions at Mn3+ ion sites where 
the charge density is reduced and there is strong FM exchange coupling between Mn^"*" ions 
lying approximately in the ac plane (-3.7 and -6.0 meV, UHF calculations Table [VIII|) . 

The three Mn-Mn interactions described by Goodenoughil are now tentatively identified 
in cubic and Jahn- Teller distorted LaMnOa using charge densities on the sites and 
orbital ordering at the neighbouring Mn ion sites. Type I interactions are found for: d^2_y2 
da,2_y2 orbital ordering in cubic LaMnOa for both Jy = 5.1 meV and = 7.2 meV (Fig. |^, 
middle panels); d^^^^r^ d^z^-r^ orbital ordering for = 14.2 meV (Fig. ^ top right panel). 

Type II interactions are found for: dx2-y2 d^z'^^r'^ orbital ordering in cubic LaMnOs 
(Fig. ^ bottom panels) for both and Jy = -6.0 meV; d^^'^^r'^ d3y2_r2 orbital ordering 
in Jahn- Teller distorted LaMnOs (Fig. ^, left panel (Jy = -8.1 meV). Charge densities are 
characterised by breaking of symmetry of the ion along the Mn-O-Mn axis. Obviously 
this can only occur when the orbital orderings on adjacent Mn ions differ, however this 
observation is worth making as such symmetry breaking is characteristic of FM exchange 
coupling. 

Type III interactions are found for d2,z2-r2 d2,z2_r^ orbital order in cubic LaMnOa (Fig. |^, 
top left panel) where Jy = -0.1 meV; for J^ =-0.6 meV in the Jahn- Teller distorted structure 
(Fig. ^, right panel). In both cases the weak exchange coupling is associated with strong, 
quadrupolar 0^~ charge density deformation. 

Cluster CI calculations provide detailed information on the exchange coupling mecha- 
nism. Fundamental SAF singlet states for clusters representing both CaMnOs and LaMnOa 
lie above the fundamental SAF high spin multiplicity states; this is expected to be the case 
for a wide range of magnetic ions exchange coupled via a closed shell anion. The ground state 
for the pair of magnetic ions is AF when the difference in correlation energies of the singlet 
and high spin multiplicity states exceeds the singlet/high spin state splitting, otherwise it is 
FM. Correlation energies for singlet states exceed correlation energies of the corresponding 
high spin multiplicity states in the three cases studied here (Fig. This is also likely to be 
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true for a wide range of magnetic ions which are exchange coupled via a closed shell anion 
as there are many more singlet SAP's than high spin SAP's in any particular active space. 
Por example, in the active space used for the LaMnOs cluster CI calculations there are over 
18,000 singlet SAP's compared to over 1,500 nonet SAP's, which simply reflects the fact 
that there are many more ways to arrange spin-coupled electrons to form singlet states than 
there are to form nonet states for a specific number of electrons. Only a few of either the 
singlet or nonet SAP's appear in the ground state wave functions with a signicant weight, 
but since there are so many more singlet than nonet SAP's, it is not surprising that the 
singlet state correlation energy is larger. 

When one empty Mn e^ orbital is available to accept one or two electrons from an 
ion, as is the case for Jy in Pnma LaMnOs, the singlet state correlation energy is only 
slightly larger than the nonet state correlation energy (82.8 versus 79.7 meV) and the nonet 
state is the ground state. However, when two empty Mn e^ orbitals are available, as is the 
case for in Pnma LaMnOs and J in CaMnOa (Pig. |^), singlet state correlation energies 
are significantly larger than the nonet(septet) state correlation energies (105.4 versus 83.3 
meV (LaMnOs) and 153.2 versus 133.4 meV (CaMnOs)) and the ground states are singlets. 

Exchange coupling in CaMnOs and LaMnOs is largely due to quantum fluctuations in 
the ground state in which one electron is excited from an 0^~ ion into an Mn e^ orbital. 
Pluctuations in which an electron is exchanged between t2g orbitals enter the singlet state in 
MnaOn" and MusOi?" clusters but are not the main contributors to the exchange interac- 
tion. Parallel studies of exchange coupling in La2CuO40 show that da;2_y2 dx2-y2 exchange 
interactions dominate the exchange coupling in La2Cu04 and O to dj;2_y2 excitations have 
a lesser weight than Cu^^Cu'^"'' excitations in the La2Cu04 ground state. This difference in 
exchange coupling mechanism most likely reflects trends in effective Hubbard U parameters 
for Mn^+Z^^ and Cu^+ and O 2p to metal d excitation energies. 
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APPENDIX A: DETAILS OF CALCULATIONS 

The methods used to generate localised orbitals and point charge arrays and the basis 
sets and CI computer codes used in this work are described in this section. 

UHP calculations on the crystalline solid were performed using the crystal98 code0. 
The basis sets used for both crystal UHP calculations and cluster CI calculations were 
identical Gaussian orbital basis sets designed for use in the solid state. They are slightly 
modified versions of the basis sets available from the crystal98 websiteB. Outer exponents 
of the Gaussian functions were modified so that the total energy in a UHP calculation on 
CaMnOs was minimised. The original basis sets had been optimised for different Mn ionic 
solids. The basis sets used in all calculations are the all electron basis sets for: Mn (86- 
411d41Gii with two d orbital exponents, optimised for CaMnOs by changing the outer d 
exponent to 0.259); O (8-411G@ with principal quantum number up to n = 4, optimised 
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for CaMnOs by changing the outer sp exponents to 0.4763 and 0.22); Ca (86- 511d3G@ with 
the outer d orbital exponents optimised for CaMnOs to 3.191, 0.8683 and 0.3191) and an 
La basis set optimised for the La^^ ion@. The La basis set used in this work differs from the 
cited basis in that the 5d orbital was removed from the basis and the 6sp and 7sp orbitals 
were replaced by a single sp orbital exponent of 0.3917. 

Different orbital ordered states in UHF calculations were obtained using a feature in the 
CRYSTAL98 code which increases the diagonal element of the Fock matrix corresponding to 
a particular orbital for the first few iterations of the calculation. This results in that orbital 
being unoccupied during those SCF cycles and allows the wave function to converge to a 
state which is a local energy minimum with a particular orbital ordering. 

High spin multiplicity states, such as the septet and nonet states of the clusters used here, 
are generally well described by a self consistent field (SCF) restricted open shell Hartree-Fock 
(ROHF) wave function. All electrons on the cluster were treated explicitly - no pseudopo- 
tential approximation was used, except in the test calculation with a La^+ pseudopotential 
described above. CI calculations were performed in localised orbital bases, rather than the 
canonical molecular orbital bases derived from the SCF ROHF calculations. Localisation 
of SCF ROHF molecular orbitals was performed using the Foster-Boys algorithm^, which 
generates localised orbitals with maximally separated centroids. Doubly occupied O 2p or- 
bitals, singly occupied Mn t2g (or t2g and an e^ orbital for Mn20ii^) orbitals and unoccupied 
Mn eg orbitals were localised in three separate localisation steps. These must be performed 
separately in order to preserve invariance of the ROHF total energy, since any mixing be- 
tween orbitals of different occupancy will increase the total energy. In the localised orbital 
ROHF wave functions for either spin state of the Mn20}^~ and Mn20j5~ clusters, each Mn 
d electron occupies a separate orbital. 

Calculations on low spin multiplicity states of the clusters used the same sets of localised 
orbitals. They demonstrate that the localised orbitals generated for the high spin multi- 
plicity states are very good approximations to the optimal orbitals for open shell low spin 
multiplicity states and that a high spin multiplicity ROHF wave function ought to be an 
excellent starting point for perturbative calculations on high and low spin multiplicity states 
in the solid state. 

In a CI calculation on a cluster of this size it is essential to partition the orbital space into 
a core space (with doubly occupied orbitals), an active space of orbitals which are occupied or 
unoccupied in the ROHF main configuration and a space of redundant, unoccupied orbitals 
which are not used in the calculation. The active orbitals in this work were the three O 
2p localised orbitals on the central 0^~ and a set of t2g and e^ orbitals on each Mn ion. 
The Direct Multi-Reference CI mo duidil m the GAMESSIH programme was used for this 
work. The active space consisted of either the (single) fundamental SAF orbitals or the 
fundamental SAF plus all possible single or double excitations which can be made from the 
fundamental SAF into empty active space orbitals. 

Calculations were performed for clusters with no surrounding point charges and with 
point charges in a spherical volume surrounding the cluster. The radius of the sphere was 
over 20 A and included around 3300 charges. The charges were located on the ionic sites 
of either CaMnOa or LaMnOa. Mulliken populations derived from UHF crystal calculations 
were used as guides for point charge magnitudes. For CaMnOs, UHF Mulliken populations 
were Ca+^-^^Mn+^-^^Os However, in the SCF ROHF cluster calculation, this choice 
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of point charges results in Mulliken populations of Mn+^-^°0 i-3i]y[n+2.60 ^^ic central 
Mn-O-Mn chain in the cluster. The Mn and Ca point charge magnitudes were adjusted to 
Ca^-^^Mn2-S''0-^-33 and this resulted in Mulliken populations of Mn+^-^^Q-^-^^Mn+^-i^ on the 
central Mn-O-Mn chain and populations of -1.64 and -1.67 on the other two O types in the 
cluster. Note that this adjustment leaves each point charge unit cell charge neutral and the 
point charge sphere radius is adjusted so that the entire cluster has a charge near zero. The 
major changes which occur on adjusting the point charges are: charge is transferred from the 
outer O ions in the cluster to the Mn ions and central O ion, each gaining about 0.4e; the AF 
exchange constant changes from 21.0 meV to 8.1 meV, in agreement with other calculation 
methods and in reasonable agreement with experiment; the degree of correlation in the 
wavefunction decreases sharply. When a Mn20|f~ cluster with no external point charges 
is used, the Mulliken populations on the central Mn-O-Mn chain are Mn"'"^'^^0~°'^^Mn"'"^'^^ 
and the exchange coupling energy is 57 meV. 

A similar adjustment of point charge magnitudes was used for the Mn20ii~ cluster 
CI calculations. The MuUiken populations determined from UHF crystal calculations on 
the experimental Pnma structure were La"'"^'^^0~^'^^'~^'^^Mn"'"^'^'*. Cluster point charges of 
La+2-^°0-i-^°'-^-s°Mn+2-6° resulted in Mulliken populations of Mn+^-^^O"^-^^ in the Pnma 
structure. 
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FIGURES 

FIG. 1. Pnma structure of LaMnOs according to ElemansS. Mn-0 bonds are shown explicitly. 
Mn ions are dark spheres, O ions are light spheres and La ions are unconnected light spheres. Mn 
ions labelled 1 and 2 are AF coupled (J_l) and Mn ions labelled 2 and 3 are FM coupled (J||). The 
cluster used to compute the AF coupling constant had the same structure as Mn ions 1 and 2 and 
their associated ion quasi-octahedra. The cluster used to compute the FM coupling constant 
had the same structure as Mn ions 2 and 3 and their associated O^^ ion quasi-octahedra. 

FIG. 2. Empty orbital ordering in LaMnOs. (A) The empty orbital arrangement which results 
when occupied orbitals are dr^2_y2 dx2_y2 ordered. (B) Empty orbital arrangement with AF spin 
coupling favoured, (C) empty orbital arrangement with FM spin coupling favoured, (D) empty 
orbital arrangement with weak spin coupling. 

FIG. 3. Localised orbital basis used for CaMnOs cluster CI calculations. Top panel: O 2p^ 
orbital; middle panel: Mn dxz orbital; bottom panehMn d2z^_j.2 orbital. The latter is the empty 
eg orbital responsible for exchange coupling. 

FIG. 4. Localised orbital basis used for exchange constant cluster CI calculation for 
LaMnOs. Top panel O 2pz orbital; middle panel: filled e^ orbital perpendicular to Mn-O-Mn 
axis; bottom panel: empty e^ orbital oriented along Mn-O-Mn axis. 

FIG. 5. Exchange coupling constants for CaMnOs and LaMnOs from CI cluster calculations 
with varying Mn ion Mulliken populations. The variation in Mulliken population was induced by 
changing the magnitude of point charges at Mn and La or Ca ion sites. 

FIG. 6. Correlation energies in CaMnOs and LaMnOs. Magnitudes of correlation energies are 
illustrated by vertical arrows and are given in meV in plain text. Energy differences between low 
and high multiplicity spin states are given in italics and occupancies of the fundamental SAF in 
each state are given at the base of each arrow. The horizontal line is the SCF ROHF energy for 
each state. 

FIG. 7. Charge density difference plots for cubic LaMnOs with (a,b) d^z^_j.2 ds22_^2 , (c,d) 
da,2_y2 d^.2_y2 , (e,f) dx2_y2 ds22_r2 orbital ordering. Left and right panels show density differences 
in the xy and xz planes, respectively. The differences in charge densities are the UHF SCF density 
for the solid minus the UHF SCF densites for the 0^~ ions and the Mn^+ ion. 

FIG. 8. Total energies of d^2_y2 dx2_y2 , ds^2_^2 ds22_,.2 and dx2_y2 d3^2„j,2 orbital ordered 
structures with G-AF magnetic ordering. The reference energy is the cubic LaMnOs da,2_y2 d^2_y2 
FM energy (Table 0). 
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FIG. 9. Charge density difference plots for LaMnOs with a 5 percent Jahn-Teher distortion in 
the xy plane. Panels (a) and (b) show density differences in the xy and xz planes, respectively. 
The differences in charge densities are the UHF SCF density for the solid minus the UHF SCF 
densites for the ions and the Mn^+ ion. 
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TABLES 

TABLE I. Relative energy and magnetic moment per Mn ion in CaMnOs. 



Spin Ordering Relative Energy (meV)*^ /^(/^b) 

FM 00 sm 

A-AF -23.7 3.27 

C-AF -45.1 3.15 

G-AF -64.3 3.23 



^Lattice constant 3.73 A 

TABLE II. Exchange constants in CaMnOa derived from experiment and ab initio calculations. 

J(meV) 

Experiment*^ 6.6 
Cluster CI^ 8.1 
Model Hamiltonian'^ 9.5 
UHF^ 10.7 



'^Rushbrooke et al.a 

^This work. Lattice constant 3.73 A 

'^Meskine et a/.i 



TABLE III. Relative energy and SAF occupation numbers for singlet and septet states of 
Mn20|^^ cluster representing CaMnOs. 

State Energy (meV)'' Main SAF t2g Exchange O to e^ (le) O to e^ (2e) 

singlet +3^6 1.0000 0.0000 0.0000 0.0000 

septet'" 1.0000 0.0000 0.0000 0.0000 

singlet<= -149.6 0.9926 0.0005 0.0038 0.0017 

septet^ -133.4 0.9943 0.0000 0.0027 0.0017 

^Energies are relative to the restricted open shell Hartree-Fock septet state 
'"Fundamental SAF only 

'^Fundamental SAF + all single and double excitations in active space from fundamental 
SAF 
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TABLE IV. Structural parameters in Jahn- Teller distorted LaMnOa and Pnma LaMnOa de- 
termined by experiment and total energy minimisation. Each cell is a \/2x2x\/2 doubling of the 
primitive perovskite unit cell. 



Ion 


X 


y 


z 


La*^ 


0.549 


0.250 


0.010 


Mn 


0.000 


0.000 


0.000 





-0.014 


0.250 


-0.070 





0.309 


0.039 


0.224 


La^ 


0.517 


0.250 


0.001 


Mn 


0.000 


0.000 


0.000 





-0.002 


0.250 


-0.027 





0.290 


0.014 


0.237 


La'^ 


0.500 


0.250 


0.000 


Mn 


0.000 


0.000 


0.000 





0.000 


0.250 


0.000 





0.2625 


0.000 


0.2625 



^Experimental structurdll (Fig. |I|), lattice parameters a=5.742A, b=7.668A, c=5.532A 
^Optimised structure, lattice parameters a=5.740A, b=7.754A, c=5.620A 



Jahn- Teller distorted structure, lattice parameters a=5.590A, b=7.905A, c=5.590A. Note 
that the Jahn- Teller distortion is in the xz plane in this Table to allow easy comparison 
between its structural parameters and those of the Pnma structures. Elsewhere in this work 
the Jahn- Teller distortion is assumed to be in the xy plane. 



TABLE V. Relative energy and magnetic moment per Mn ion in cubic LaMnOs with various 
spin and orbital orderings. 

Spin and Orbital Ordering*^ Relative Energy (meV) /^(/^b) 



FM (lri.2_y2 da;2_y2 


0.0 


4.00 


FM d^2_y2 d3^2_j,2 


-131.5 


4.00 


FM d3^2_j,2 d322_^2 


-6.1 


3.99 


A-AF d,j.2_y2 6.^2 _y2 


-14.4 


4.05 


A-AF d^z2_j.2 d322_^2 


-34.4 


3.87 


G-AF d2,.2_j^2 6.^2 _y2 


-34.9 


3.88 


G-AF d^2_y2 d^^2_^2 


-95.4 


3.89 


G-AF d^z2_j.2 d3^2_j.2 


-34.0 


3.88 



^Lattice constant 3.934A 
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TABLE VI. Exchange constants in cubic LaMnOa with various orbital orderings. 



Spin and Orbital Ordering^ 


(meV)^ 


J|| (meV)^ 


d-x2—y2 djp._y2 


7.2 


5.1 


dx^—y^ d322_^2 


-6.0 


-6.0 


d3^2_r2 d^z2_r2 


14.2 


-0.1 



'^Lattice constant 3.953A 

^Exchange constant for Mn ions coupled perpendicular to ac plane 
*^Excliange constant for Mn ions coupled parallel to ac plane 



TABLE VIL Relative energy and magnetic moment per Mn ion in Pnma and Jahn- Teller 



distorted cubic LaMnOa. 



Structure and Spin Ordering 


Energy(meV) 




Pnma (Experiment) FM'^ 


0.0 


4.00 


Pnma (Experiment) A-AF 


-1.2 


4.00 


Pnma (Experiment) G-AF 


13.9 


3.96 


Pnma (Optimised) FM*^ 


0.0 


4.00 


Pnma (Optimised) A-AF 


-2.0 


3.96 


Pnma (Optimised) G-AF 


21.9 


3.94 


Jahn-Teller FM= 


0.0 


4.00 


Jahn-Teller A-AF 


1.1 


3.98 


Jahn-Teller G-AF 


33.6 





^Reference energy is 194 meV above optimised FM Pnma structure (Table ^V]) . 
^Reference energy is that of this structure and magnetic order (Table |IV|). 
'^Reference energy is 8 meV below optimised FM Pnm,a structure (Table ^Vf ). 
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TABLE VIII. Exchange constants in Pnma LaMnOa derived from experiment and ab initio 
and model Hamiltonian calculations. 





J_L (meV)*^ 


J|| (meV)^ 


Experiments^ 


4.8 


-6.7 


UHF (Experiment ) 


0.6 


-3.7 


UHF(Experiment)'^ 


0.8 


-3.5 


UHF(Optimised)^ 


1.0 


-6.0 


UHF(Jahn-Teller)g 


-0.6 


-8.1 


LSDA(Experiment)'^ 


-3.1 


-9.1 


Cluster CI(Experiment)'^ 


3.3 


-3.6 


Cluster Cl(Optimised)'' 


5.1 


-7.4 


CI(Optimised/La pseudopotential)^ 


5.2 


-7.4 


model Hamiltonian' 


4.8 


-3.7 



'^Exchange constant for Mn ions coupled perpendicular to ac plane 

'^Exchange constant for Mn ions coupled parallel to ac plane 

'^Hirota et a/.i;Moussa et a/.i 

'^This work. Elemans structurelii(Table ^Vl) - 

^Su et alE. 

^This work. Optimised structure (Table |^). 

^This work. Jalin- Teller distorted structure (Table fV\ j. 

^Soloyev et alB. 

'Meskine et alM 
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TABLE IX. Relative energy and orbital occupation numbers for singlet and nonet states of 



MnaOl?" 


cluster representing 

IT C 


; LaMnOs. 








State 


Energy (meV)'^ 


Mam SAF 


t2g Exchange 


(J to Cy (lej 


to Cg (2e) 


smglet 


+11.9 


1.0000 


0.0000 


0.0000 


0.0000 


nonet 


0.0 


1.0000 


0.0000 


0.0000 


0.0000 


singlet'^'^ 


-93.5 


0.9937 


0.0006 


0.0037 


0.0007 


nonet'^'^ 


-83.3 


0.9954 


0.0000 


0.0030 


0.0007 


singlet*^ *^ 


+17.9 


1.0000 


0.0000 


0.0000 


0.0000 


nonet'^*^ 


0.0 


1.0000 


0.0000 


0.0000 


0.0000 


singlef^^ 


-64.9 


0.9949 


0.0004 


0.0025 


0.0006 


nonef^^ 


-79.9 


0.9946 


0.0000 


0.0038 


0.0008 



^Energies are relative to the restricted open shell Hartree-Fock septet for the whole cluster, 
i.e. per Mn ion pair. 

'^Fundamental SAF only 

'^Fundamental SAF + all single and double excitations in active space from fundamental 
SAF 

'^Jj, calculation 
^Jii calculation 
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Mn ion charge 



CaMn03 



septet singlet 



LaMn03 
interplanar intraplanar 
nonet singlet nonet singlet 



133.4 



16.2 

0.9942 ^ 

0.9926 



153.2 



83.3 

70.2 



105.4 



79.7 



82.8 



14.7 
0.9975 



0.9977 



0.9973 



0.9968 



